Weak Alfven- Wave Turbulence Revisited 
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Weak Alfvenic turbulence in a periodic domain is considered as a mixed state of Alfven waves 
interacting with the two-dimensional (2D) condensate. Unlike in standard treatments, no spectral 
continuity between the two is assumed and indeed none is found. If the 2D modes are not directly 
forced, k~ 2 and fc" 1 spectra are found for the Alfven waves and the 2D modes, respectively, with the 
latter less energetic than the former. The wave number at which their energies become comparable 
marks the transition to strong turbulence. For imbalanced energy injection, the spectra are similar 
and the Elsasser ratio scales as the ratio of the energy fluxes in the counterpropagting Alfven waves. 
If the 2D modes are forced, a 2D inverse cascade dominates the dynamics at the largest scales, but 
at small enough scales, the same weak and then strong regimes as described above are achieved. 
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Introduction. It has been understood for many years 
that small-scale turbulence of a conducting fluid or 
plasma in a strong magnetic field consists of Alfven- wave 
packets This is true in most astrophysical plasmas, 
including the weakly collisional ones, where Alfvenic fluc- 
tuations populate the scales above the ion Larmor scale 
Q. On theoretical 0-(l], numerical and observa- 

tional grounds, it appears clear that Alfvenic tur- 

bulence is anisotropic with k± 3> kn . The parallel scales 
are associated with the propagation of Alfven waves, the 
perpendicular ones with the nonlinear interaction be- 
tween them. The relative importance of these effects de- 
pends on the corresponding time scales, ta ~ (^|| u a) _1 
and T n i ~ (k±u±)~ 1 , where va is the Alfven speed and uj_ 
the perpendicular velocity perturbation. When ta S> r n i, 
the nonlinearity dominates and the turbulence is effec- 
tively two-dimensional (2D); when ta <C t u \, the wave 
propagation dominates and the turbulence is weak. 

A causality argument suggests that a pure 2D regime 
cannot be sustained: for any given kj_, motions in two 
planes perpendicular to the mean field and separated by 
a distance ~ fc^ 1 can only remain correlated if the time 
it takes an Alfven wave to propagate between the planes 
is ta < T n \. Thus, an initially 2D perturbation will nat- 
urally decay into a state of "critical balance," ta ~ t u \. 
It has been argued , |T3l that weak turbulence, the 
limit case opposite to 2D, will also approach critical bal- 
ance via a perpendicular cascade in which r n i becomes 
ever smaller until r n \ ~ ta at some sufficiently small 
scale. The critical balance thus appears to be the fun- 
damental physical principle underpinning Alfvenic tur- 
bulence (and possibly in general turbulence in systems 
that support propagation of waves [H|). However, the 
structure of critically balanced turbulence remains con- 
tentious and poorly understood. Efforts to improve this 
understanding have often turned to various insights from 
the theory of weak turbulence for intuition and guid- 



ance. Weak Alfven-wave turbulence itself, while having 
for some time enjoyed the reputation of a solved problem 
[l6l | , has nevertheless recently been subject of several new 
investigations that amended or disagreed with the estab- 
lished paradigm [I9I - I22I ] . This paper is a contribution 
to this revisionist tendency, focusing on the structure of 
weak Alfven-wave turbulence in finite periodic domains 
and on the transition to critical balance. While a peri- 
odic box may be an artificial setting, it is ubiquitous in 
numerical experiments. It is, therefore, important to un- 
derstand turbulence in such domains and the extent to 
which it might belong to the same universality class as 
turbulence in natural systems. 

While it might appear that weak turbulence is an ana- 
lytically tractable and, therefore, easily understood limit, 
the weak turbulence of Alfven waves is, in fact, diffi- 
cult to treat in a rigorous fashion because of a special 
role played by the modes with fc|| = 0. Since Alfven 
waves have frequencies lo^ = ±k^VA and only counter- 
propagating waves can interact, the resonance conditions 
and fci + &2 = &3 imply that at least one 



mode in any interacting triad must have fcii = [23|, 
These modes are not Alfven waves but rather 2D mo- 
tions for which ta = 00, so they cannot be treated by 
the weak-turbulence approximation. The standard rem- 
edy for this complication has been to proceed with the 
weak-turbulence expansion anyway, assuming formally 
that the k\\ spectrum of the Alfven waves is continuous 
across fen = 1|| 17 1. The resulting theory predicts a 



kj_ scaling of the energy spectrum. As this scaling is cor 



roborated by direct numerical simulations [251 ] . the con 



tinuity assumption might seem to be vindicated. In this 
paper, we propose a different way of treating the fcii = 
modes, with no assumption of their spectral continuity 
with the Alfven waves. It leads to a new phenomenolog- 
ical theory of the weak Alfven-wave turbulence and to 
distinct scaling predictions for the energy spectra of the 
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Alfven waves and of the feii — modes. We discuss vari- 
ous regimes: balanced and imbalanced, containing forced 
2D motions or otherwise, and also describe the transition 
to critically balanced strong turbulence in a new way. 

Scaling theory. Let us start with the equations of re- 
duced magnetohydrodynamics (RMHD) [2f 



271, which 



can be shown to describe correctly the anisotropic 
Alfvenic fluctuations both in an MHD fluid [28j], and, 
above the ion Larmor scale, even in weakly collisional ki- 
netic plasmas In RMHD, the velocity and magnetic 
field perturbations perpendicular to the mean field Bq = 
vaz (the Alfvenic perturbations) are two-dimcnsionally 
solenoidal, so they can be expressed in terms of stream 
and flux functions: u± = z x Vj_$ and SB± = zx Vi$. 
The RMHD equations can then be written in terms of the 
Elsasser potentials C = $ ± as follows [H 

d t t ± Tv A d z ( ± = N[(T,C ± ]+F ± , (1) 

where is the stream function of a body force repre- 
senting energy injection and the nonlinear term is 

MC^C*] = -\ vi 2 ({C + , vie} + {C, viC} 

TVi {C,C}), (2) 

where {A, B} = z ■ (V±A x V±B). 

Let us now Fourier transform in z, factor out the 
oscillating in time part of the solution, ^(x,y 7 z) = 
Sfe (k n (x,y)e ±lk tt VAt+lk u z , and write separately the evo- 
lution equations for the Alfven waves (fen ^ 0), 



KG, =W,C*j | ] + JV[cr„ J Cf]e 

+■ E N ^t,-^ eTi2k " VAt 



±l„=F*2fe||«At 



and the k\\ = modes, 
dt£ = N[£,(±] + 



E^,A 

fc,,/0 



= =Fi2/c I. da* 



(3) 



(4) 



The standard approximation of the weak-turbulence 
theory is, roughly speaking, to neglect the nonlinear 
terms in Eq. ([3]) that have oscillatory factors, so the dom- 
inant effect is the Alfven-wave "scattering" off the fc|| = 
modes (the first term on the right-hand side). This gives 
rise to a cascade of energy to small perpendicular scales 
(large fcj_), while the transfer of energy from the directly 
forced ku to other fey is small. For simplicity, let us as- 
sume that only one fey = fep is forced. If we denote 
Cn(k±) the characteristic amplitudes corresponding to 
fen = nfep and the perpendicular wave number fcj^ and 
take the interactions to be local in k±, we may estimate 



(vi^Wo^.Ci] ~ fcic T (^)[Cr t (fcx)] 2 ~ £ d 



(5) 



where e ± = ((Vj_C) ' (Vj_-F ± )} is the mean power in- 
jected by the forcing into the "+" and "— " modes. We 
first consider the balanced case: e + = e~ = e 1431 . 



Eq. iJSJ) shows that in order to make scaling predictions 
for the Alfven waves, we must know the scaling of the 
amplitudes of the fey = modes — these modes deter- 
mine the cascade rate ~ fexCo^J-) for the Alfven waves. 
The = modes arc very different from the Alfven 
waves: they are described by 2D magnetohydrodynam- 
ics with an oscillatory nonlinear source term representing 
the coupling of feii and —ku Alfven waves [Eq. (3)]. We 
shall first consider the situation in which the feii = 
modes are not forced externally, so this nonlinear source 
is the only source of energy in the ku = modes. Since 
the nonlinear source in Eq. ((4]) has an oscillatory factor, 
there is a strong cancellation effect and the amplitude of 
the fen = modes can be estimated as [3] 

Co ± (fc±)~^ l fcic 1 + (^)cr(fe±), (6) 

where uja = k\\fVA- Combining Eqs. {5} and (j6|), we find 

Cftfcj.) ~ (e^) 1/4 fcl 3/2 =>> (*±) ~ {eu A ) 1/2 kl\ 

(7) 

to(k±) ~ (—) C 1 =► E£fr±) ~ — K\ (8) 
\ujaJ ua 

where £7„ (fe±) is the one-dimensional energy spectrum, 
related to the characteristic amplitudes via fej.-E!^ (fcj_) ~ 
k±[Cn(k±)] 2 - The energy injection is balanced, so the 
the "+" and "— " spectra have the same scaling. 

The phenomenological argument presented above has 
led to a prediction of the Alfven-wave spectrum [Eq. (J7J] 
that is formally the same as the prediction of the stan- 
dard weak-turbulence theory (this is only true for the 
balanced case; see below). However, the physical origin 
of this spectrum is different and the assumption of con- 
tinuity across fcn = is certainly not satisfied: 



Tnl.O 



TA 
Tnl.l 



kj_ i 



(9) 



where ta = 



J A j 



Tnl.O ~ [kj_(o(kj_)] 1 , T nlj i 



[feiCf(fe_L)]"\ and k c = (^/e) 1 / 2 . The amplitudes of 
the fey = modes are smaller than those of the Alfven 
waves as long as ta <C T n i, i.e., as long as the weak- 
turbulence limit holds. However, their spectrum is shal- 
lower than that of the Alfven waves, so the ratio between 
the amplitudes increases with k± until they become com- 
parable at k± ~ k c . At this point the turbulence becomes 
critically balanced and is no longer weak, developing a 
£-5/3 |J |7^7jJ or p erna p S fe~ 3 / 2 [3(| spectrum. 

Higher-ku modes. In the same way that a small leak- 
age of energy from the forced modes (feii = feiif) via os- 
cillatory nonlinear couplings gives rise to a spectrum of 
feii = modes, similar leakages arise from feii = fep to 
feii = 2feiif and from there onwards to fey = 3fep, 4fc||f, 
. . . , nfep (the third term on the right-hand of Eq. |3])). 
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Note that the standard weak-turbulence treatment Q 
of the imbalanced case was insufficient to fix the individ- 
ual scalings of the "+" and "— " spectra [l7j]. Our argu- 
ment does not have this problem and is able predict the 
relationship between amplitudes and fluxes [Eq. (fH)) ]. for 
which the standard weak-turbulence theory had to make 
recourse to solving the kinetic equation, formally invalid 
for the fc|| = modes [HI]. 

The feii = modes are still small compared to the 
Alfvcnic modes: 



WEAK. 



Ef{k A 



'nl,l 
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(15) 



FIG. 1: A schematic illustration of the weak-turbulence spec- 
tra and transition to strong, critically balanced turbulence. 



These modes do not contribute to the scattering of the 
driven modes and are just induced by the latter. Analo- 
gously to Eq. ([5]), it is not hard to see that for n > 2, 

^ n -\t(k±)} n 



£(*x) 



(10) 



whence follow the spectra of these modes [using Eq. Q] 



• 3(n-l) 



E n {k A 



T n/2 



,2(»-l) 1W (3n-4)/2 



kr 



(ii) 



A A 

Similarly to Eq. (JSJ, we have E n (k±) / E\(k±) ~ 
(k±/k c ) n , so the amplitudes of the fe|| = nk\\{ modes 
with n > 2 become comparable to the amplitude of the 
forced mode (n = 1) at the same wave number as does 
the amplitude of the fey = mode. This situation is 
illustrated schematically in Fig. [TJ 

Imbalanced turbulence. If the energy injection into 
the "+" and "— " Elsasser modes is not the same, say, 
e + > s~ , the resulting turbulence is known as imbalanced 
- this is, in fact, a generic situation in the solar wind 
31 , 32[ and also in numerical simulations if one considers 
local subdomains of the simulation box 33 1 . 

Our arguments are easily adapted to this case. Eqs. 
(|5l) and ([6]) hold unchanged (note that the latter formula 
Cq"). Eqs. ([7]) and flSJ) generalize to 

1/4 



implies £o" 



( £±) 3/8 ( 



Co=(*J 



(£T) l/8 k 3/2 



+ .-U/4 



1 ( £ T)l/4^ 



(12) 



(e+e 



uj a k± 



E±{k A 



1/2 



(13) 

These imply that the cross-helicity and the Elsasser ratio 
in weak imbalanced turbulence are idependent of k±: 

Elm^l. (14) 



E+(k ± )~Ec(k x ), 



where the nonlinear times r. 
different for the "+" and 



numbers fer 1 



3/2 



ni,i [ k \C\{k±)] 1 arc now 
■" modes and so are the wave 
(e ± ) 1 / 4 (e =F )- 3 / 4 at which the weak- 



turbulence approximation breaks down. It must of course 
break down already at the smaller of the two, viz., k~ 
(we have assumed e + > e~). This opens the possibility 
of a "twighlight range" of k± between k~ and fe+ (or an- 
other threshold dependent on the intermediate scalings) 
- an extended transition from weak to strong regime, 
where the "— " modes are strongly nonlinear, while the 
"+" modes are not. It is indeed possible to construct 
such mixed theories, featuring steeper spectra for the 
"+" modes and shallower ones for the "— " modes (cf. 
[34]]). We are tempted to speculate that this might help 
explain the origin of apparently non- universal (and differ- 
ent) slopes of the "+" and "— " spectra found in simula- 
tions of strongly imbalanced Alfvenic turbulence [3HI . (3|| . 
However, the more or less arbitrary assumptions neces- 
sary to fix scalings in this regime and the consequent 
uncertainties in the outcome are so numerous that we 
prefer not to treat this subject here. The precise scal- 
ings of strong turbulence in the imbalanced regime also 
remain theoretically uncertain 

Case of hydrodynamically forced feii = modes. We 
have so far considered a special case in which the fcy = 
modes were not forced. Let us now allow comparable 
power to be injected into fey = as into feii = fep. The 
situation is now radically different because the amplitude 
of the feii = modes is no longer determined by a small 
leakage of the Alfven-wave energy via the oscillatory term 
in Eq. but by a direct forcing. If we ignore the oscilla- 
tory term altogether (assuming it averages out to lowest 
order in Ta/Thi), the feii = modes decouple and form 
an independent 2D turbulent condensate. We write the 
equations for this condensate in terms of its velocities, 
given by the stream function <!>o = (Co~+Co")/2, and mag- 
netic fields, given by the flux function VPg = (Co~ — Co)/^ : 



cWi$ + {$0, Vi$ } = {* , vi* } + ViF , 



fcii/O 



(16) 
(17) 



4 




fcn k m WEAK k c STRONG ■ 

WEAK + hydro 2D 



FIG. 2: A schematic illustration of the weak-turbulence spec- 
tra for the case of velocity-forced ku = modes. 

where we have returned to the assumption that all forcing 
is in the velocity field. Therefore, the magnetic flux func- 
tion "to is a passive scalar and its only source of energy 
is the oscillatory coupling to the Alfven waves, which has 
consequently been retained in Eq. (|17j). Since this energy 
source is small, <C <£>o and the Lorentz force can be 
neglected in Eq. (fP6|) . leaving a 2D Eulcr equation. For 
the forced Alfven waves with fen = fep, we have from from 
Eq. ([3]), neglecting the oscillatory terms and ^o, 

d t Ct = N[* , Ct) + F 1 e^ v * t . (18) 

The unforced Alfven waves with ku ^ fcp will have small 
amplitudes due to the oscillatory terms in Eq. ((3|). 

Thus, the turbulence has split into the following dis- 
tinct components: a 2D (ku = 0) forced hydrodynamical 
condensate, an ensemble of forced Alfven waves passively 
advected by this 2D condensate [Eq. (fT8l) ], small 2D mag- 
netic fluctuations feeding off the Alfven waves and also 
passively advected by the 2D hydrodynamic condensate 
[Eq. (|17|)]. and small unforced Alfven waves again pas- 
sively advected by the 2D condensate and feeding off 
the forced modes and each other (in a way similar to 
that described above for the case of no forcing of the 
ku = modes). This situation is somewhat similar to 
the "slaved" regime proposed in , where everything is 
unilaterally controlled by the 2D condensate. 

Since the condensate is hydrodynamical and 2D, one 
expects an inverse energy cascade to perpendicular scales 
larger than the forcing scale (in a finite system leading 
to energy accumulation at the system scale) (HJ. Below 
the forcing scale, the direct enstrophy cascade produces 
a well-known kinetic-energy spectrum: 

Ea(k±) ~ fc±[$o(fc±)] 2 ~ 7ofcI 3 , (19) 

where 70 = (fe^eo) 1 ^ 3 is the rate of strain (the same 
for all modes), £0 = ((Vj_$o) • (Vj_Fo)) the mean power 
injected into the fcii = modes and k±{ the perpendicular 



wave number of the forcing [46( . The spectra of the forced 
Alfven waves and the fen = magnetic fluctuations follow 
from Eq. © with (± = $ and Eq. © with $ = #0, 
respectively: 

^(fc±)~fc±[^(fc±)] 2 ~-fc:\ (20) 

7o 

M (fc ± )~fc ± [* (fc ± )] 2 ~ (—) k±, (21) 

where e± is the mean power injected into the ku = fciif 
Alfven waves (assumed balanced) . This situation is illus- 
trated in Fig. [2] The intersection wave number between 
E and E\ is ko = (ll/si) 1 ^ 2 ~ fc_Lf if £0 ~ £i- The inter- 
section between Eq and Mo is at k m = (jqCja/si) 1 ^ 2 ■ For 
fcj_ > k m , the rate of strain 70 of the hydrodynamic con- 
densate is smaller than the rate of strain associated with 
the fcy = modes induced by the oscillatory coupling 
to the Alfven waves (70^1,0 < 1) and so the turbulence 
reverts to the regime described earlier for the case of no 
forcing of the ku = modes. 

Other regimes. Possibilities proliferate if we allow im- 
balanced energy injection and/or let the ku = modes 
be forced magnetically as well as hydrodynamically. The 
majority of the resulting regimes arc probably not physi- 
cally realizable, so we will not pursue this line of inquiry 
further, except for the following observation. If is 
forced, this direct injection of energy into 2D magnetic 
fluctuations will in general trump the oscillatory coupling 
to Alfven waves (the right-hand side of Eq. (fT7|) ). Thus, 
the 2D condensate decouples fully from the latter. Mag- 
netic field of the ku = modes is now dynamically signifi- 
cant, so there will be a 2D inverse cascade of the variance 
of the magnetic flux function, l^ol 2 ! and a direct cascade 
of the two Elsasser energies iV^Co^l 2 - The inverse cas- 
cade of l^ol 2 to fcj^ < k±{ produces a fc^ 1 ^ 3 spectrum 
of magnetic energy and possibly a shallower spectrum 
of the kinetic energy [39j]. This may well be the phys- 
ical mechanism responsible for the relative preponder- 
ance of magnetic energy over kinetic at large scales in 
weak MHD turbulence [2~lj |. The direct cascade is nu- 
merically known to give rise to fc^ 3 ^ 2 spectra [40l - [4^ | [47T] . 
The Alfven waves (fcii = ku{) will be passively mixed by 
this strong 2D MHD turbulence. However, by the causal- 
ity argument given in the Introduction, the latter is, in 
fact, likely to be unstable, develop parallel decorrelations 
and become strong, 3D and critically balanced. 

This work was supported by STFC (AAS and TAY) 
and the Lcvcrhulmc Trust Network for Magnetised 
Plasma Turbulence. Much of it was completed during 
the Programme "Frontiers in Dynamo Theory" at Insti- 
tut H. Poincare, Paris (March- April 2009), whose hospi- 
tality AAS and SVN gratefully acknowledge. 



5 



* Electronic address: a.schckochihml@physics.ox.ac.uk 
[1] R. H. Kraichnan, Phys. Fluids 8, 1385 (1968). 
[2] A. A. Schekochihin et al., Astrophys. J. Suppl. 182, 310 

(2009) . 

[3] D. Montgomery and L. Turner, Phys. Fluids 24, 825 
(1981). 

[4] P. Goldreich and S. Sridhar, Astrophys. J. 438, 763 

(1995) . 

[5] P. Goldreich and S. Sridhar, Astrophys. J. 485, 680 
(1997). 

[6] J. V. Shebalin et al., J. Plasma Phys. 29, 525 (1983). 
[7] J. Cho and E. T. Vishniac, Astrophys. J. 539, 273 (2000). 
[8] J. Maron and P. Goldreich, Astrophys. J. 554, 1175 
(2001). 

[9] W.-C. Miiller et al., Phys. Rev. E 67, 066302 (2003). 
[10] J. W. Bieber et al., J. Geophys. Res. 101, 2511 (1996). 
[11] A. Bigazzi et al., Astrophys. J. 638, 499 (2006). 
[12] L. Sorriso-Valvo et al., Europhys. Lett. 75, 832 (2006). 
[13] T. S. Horbury et al., Phys. Rev. Lett. 101, 175005 (2008). 
[14] R. T. Wicks et al., Mon. Not. R. Astron. Soc. 407, L31 

(2010) . 

[15] C. H. K. Chen et al., Mon. Not. R. Astron. Soc. 415, 
3219 (2011). 

[16] S. Galtier et al., J. Plasma Phys. 63, 447 (2000). 
[17] Y. Lithwick and P. Goldreich, Astrophys. J. 582, 1220 
(2003). 

[18] S. V. Nazarenko and A. A. Schekochihin, J. Fluid Mech. 

677, 134 (2011). 
[19] S. Nazarenko, New J. Phys. 9, 307 (2007). 
[20] S. Boldyrev and J. C. Perez, Phys. Rev. Lett. 103, 

225001 (2009). 

[21] Y. Wang et al., e-pr int larXTv:1106.2238l (2011). 

[22] A. Alexakis, e-print larXiv : 1109.2039 (2011). 

[23] D. Montgomery and W. Matthaeus, Astrophys. J. 447, 

706 (1995). 

[24] C. S. Ng and A. Bhattacharjee, Astrophys. J. 465, 845 

(1996) . 

[25] J. C. Perez and S. Boldyrev, Astrophys. J. 672, L61 
(2008). 

[26] B. B. Kadomtsev and O. P. Pogutse, Sov. Phys. JETP 

38, 283 (1974). 
[27] H. R. Strauss, Phys. Fluids 19, 134 (1976). 
[28] D. C. Montgomery, Phys. Scripta T2/1 (1982). 
[29] A. Beresnyak, Phys. Rev. Lett. 106, 075001 (2011). 
[30] S. Boldyrev, Phys. Rev. Lett. 96, 115002 (2006). 



[31] C. Tu et al., Geophys. Res. Lett. 17, 283 (1990). 
[32] E. A. Lucek and A. Balogh, Astrophys. J. 507, 984 
(1998). 

[33] J. C. Perez and S. Boldyrev, Phys. Rev. Lett. 102, 

025003 (2009). 
[34] B. D. G. Chandran, Astrophys. J. 685, 646 (2008). 
[35] A. Beresnyak and A. Lazarian, Astrophys. J. 682, 1070 

(2008). 

[36] A. Mallet et al., in preparation (2011). 
[37] Y. Lithwick et al., Astrophys. J. 655, 269 (2007). 
[38] J. J. Podesta and A. Bhattacharjee, Astrophys. J. 718, 
1151 (2010). 

[39] D. Biskamp and U. Bremer, Phys. Rev. Lett. 72, 3819 
(1993). 

[40] D. Biskamp and H. Welter, Phys. Fluids B 1, 1964 (1989). 
[41] D. Biskamp and E. Schwarz, Phys. Plasmas 8, 3282 
(2001). 

[42] C. S. Ng et al., Phys. Plasmas 10, 1954 (2003). 

[43] Physically speaking, this is the only case we are allowed 
to model by Eq. JTJ because we must in fact set F + = F~ 
in order to ensure that only velocity field is directly 
forced. Indeed, if F + - F~ ^ 0, we would be introducing 
an inhomogeneous forcing term in the induction equa- 
tion, leading to unphysical breaking of the magnetic flux 
conservation at the forcing scale. This said, F + ^ F~ 
is a popular modeling choice for simulating imbalanced 
MHD turbulence and we will consider this case later on. 

[44] A reader who does not find this estimate obvious on 
dimensional grounds, may be convinced by the follow- 
ing argument. Consider a model stochastic equation for 
the fc|| = modes: d t (o = -(o/r c + x{t)e' ibJAt , where 
t~c ~ Tni.o is the decorrelation time due to nonlinear 
mixing of £o by and the % term models the oscillatory 
coupling to the Alfven waves [see Eq. Q]. The forcing 
amplitude x(^) ~ ^ICi^Cr a ls° nas decorrelation time 
t c because ^ are also mixed by £o [see Eq. ©]. Tak- 
ing (x(*)x* "(O) = (lx| 2 ) e_ '' _ * ^ Tc an d solving our model 
equation, we obtain (|Co| 2 ) ~ (|x| 2 )/ w2 t> q.e.d. 

[45] An inverse cascade in weak MHD turbulence has re- 
cently been reported [53]. The regime with an inverse 
cascade and spectra similar to Fig. [5] was also found by 
T. A. Yousef in unpublished numerical work (2008). 

[46] In a finite box, the energy will accumulate at the box 
wave number fcxbox and so 70 = fc±box(eoi) 1//2 - 

[47] A possible explanation for this in the form of an adapta- 
tion of the dynamical alignment argument [3(| to 2D is 
found at the end of Sec. 4 of fisl ]. 



